"Single-cycle" ionization effects in laser-matter interaction 
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Abstract 

We investigate numerically effects related to "single-cycle" ionization of dense 
matter by an ultra-short laser pulse. The strongly non-adiabatic response of 
electrons leads to generation of a MG steady magnetic field in laser-solid 
interaction. By using two-beam interference, it is possible to create periodic 
density structures able to trap light and to generate relativistic ionization 
fronts. 

I. INTRODUCTION 

In the adiabatic field ionization regime, the ionization rate grows sharply when the 
electric field approaches the barrier suppression (BS) limit, i.e. when the laser intensity is 
high enough that the electron in the ground state is able to "classically" escape the atomic 
potential barrier. The ionization rate for such field strength may become higher than the 
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laser frequency and a regime in which most of the ionization is produced within a single 
laser half-cycle is achievable. 

Here we present a numerical study of some effects of ultrafast ionization in the interaction 
of a short laser pulse with an initially transparent dense medium. First, we will discuss the 
generation of megagauss steady magnetic fields in the surface "skin" layer of "solid" targets, 
i.e. slabs of hydrogen atoms with a number density close to that of a solid medium (Macchi 
et al. 1999). 

Second, we will describe effects related to the combination of two-beam interference with 
ultrafast ionization. We will show how it is possible to take advantage of this feature to 
create a layered dielectric-conductor structure able to trap the electric field, as well as a 
relativistic ionization front (Conejero Jarque et al. 1999). 



The generation of steady currents and magnetic fields by ultrafast ionization is due to 
the non-adiabatic nature of the response of initially bound electrons to a strongly ramping 
laser field. Using a following "simple-man's" model (SMM), very similar to the SMM used 
in studies of above-threshold ionization and harmonic generation in atoms, it can be shown 
that a single electron subject to an external sinusoidal intense field can acquire a steady 
velocity (Macchi et al. 1999) 



where vj is the ejection velocity of the electron, v qo = eEyo/mw, E yo is the maximum field 
amplitude and Et is the field amplitude at the instant of ionization, which will be close to the 
threshold field for barrier suppression (for hydrogen, E T ~ 0.146E aM = 7.45 x 10 8 V cm -1 , 
being E au = 5.1 x 10 9 V cm -1 the atomic field unit). 

If most of the electrons in the medium are ionized at the same instant, as may happen 
with a pulse which sharply rises above E T , one gets a net steady current which in turn 



II. GENERATION OF STEADY MAGNETIC FIELDS 
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generates a magnetic field. To obtain a larger current one may think to "tune" appropriately 
Et and E y0 . This is possible if the ionization is no longer correlated with the oscillating 
field, i.e., it is produced independently of the field itself, like in the case studied by Wilks 
et al. (1988), in which a steady magnetic field B st ~ E y o can be obtained in a very dense 
medium. For intense lasers (/ > 10 18 W cm -2 ), such a magnetic field would get values 
exceeding 100 MG and could explain (Teychenne et al. 1998) the experimental observation 
of high transparency of thin foil solid targets to 30 fs, 3 x 10 18 W cm -2 pulses (Giulietti et 
al., 1997). However, it is questionable whether this high magnetic field may be obtained with 
superintense laser pulses. In this case, in fact, the "source pulse" itself ionizes the medium 
and thus this will impose a constraint on the phase mismatch between the field and the 
velocity of the electrons. We will show by numerical simulations that the steady magnetic 
field exists but has values around 1 MG, being therefore too weak to allow enhanced laser 
propagation. 

A. PIC simulations 

First we review the results of 1D3V PIC simulations with field ionization included. We 
choose pulses with a "sin 2 " envelope and with a "square" envelope. For all the PIC runs, 
the laser frequency was uj^ = 2 x 10 15 s _1 , close to that of Nd and Ti:Sapphire lasers. The 
thickness of the target was 0.09 /zm and the density was n Q = 6.7 x 10 22 cm" 3 {oj po ~ 7loi). 
For the ionization rate we used a semi-empirical formula obtained from atomic physics 
calculations (Bauer and Mulser, 1999). The laser energy loss due to ionization is included 
introducing a phenomenological "polarization" current (Rae and Burnett 1992, Cornolti et 
al. 1998, Mulser et al. 1998). 

Fig.^ shows the spatial profiles of the magnetic field and the free electron density five 
cycles after the end of a five cycles long (At^ = 15 fs) pulse, for three different field intensities 
in the "sin 2 " shape case, and for the square profile case at the intermediate intensity value. 
The steady field is generated at the beginning of the interaction and is always much weaker 
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than the laser field, even for the most intense case (corresponding to an intensity of 3.5 x 
10 18 W cm -2 ); its sign varies according to the phase of the laser half cycle where most of the 
ionization occurs. The ionization at the left boundary is nearly instantaneous; however, even 
if the target is only 0.1 A thick, it is not ionized over its whole thickness due to instantaneous 
screening, except for the maximum intensity case. 

The fact that the produced magnetic field is much less than expected may be attributed 
to the instantaneous screening of the EM wave due to the ultrafast ionization. In fact, it is 
too weak to affect self-consistently the refractive index and as a consequence it cannot lead 
to magnetically induced transparency as hypothesized by Teychenne et al. (1998). 

B. Boltzmann simulations 

To yield a further insight into the magnetic field generated by ultrafast ionization we 
look at the results of ID and 2D Boltzmann simulations. This corresponds to the "direct" 
numerical solution of the Boltzmann equation for the electron distribution function f e = 
/ e (x, v, t), over a phase space grid: 

eE 

dtfe + V • V/ e «9 V / e = ^(i?K(x, t) 9 (v, E(x, *)). (2) 

m 

Here n a is the density of neutral atoms (supposed at rest for simplicity) and u I is the 
ionization rate. The term g(v; E) gives the "instantaneous" distribution of the just ionized 
electrons, which is supposed to be known from atomic physics. A semiclassical picture 
which allows to define and evaluate g(v; E) was given by Cornolti et al. (1998). With 
respect to PIC simulations, the Boltzmann approach has the disadvantage of larger memory 
requirements, but the advantages of reduced numerical noise and the possibility to take into 
account the full kinetic distribution of the ionized electrons. 

We first look at 2D2V Boltzmann simulations. We take a 0.25 jum, 10 16 W cm -2 laser 
pulse impinging on a solid hydrogen target with number density 2 x 10 23 cm" 3 = 12.5n c , and 
thickness 0.1 /im. The time envelope of the laser pulse is Gaussian with a FWHM duration 
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of 2 cycles. The laser spot is also Gaussian with a FWHM of 2 //m. Fig.|2| (a) shows the 
magnetic field and the density contours after the end of the laser pulse. The steady magnetic 
field has constant (negative) sign over its extension. Its maximum intensity is about 3 MG. 
Fig.|2] (b) shows the electron current density j y at the same time of the right plot of fig.0 (a). 

Among the parameters of our simulations, the magnetic field appears to be most sensitive 
to the temporal profile of the laser pulse, achieving its maximum value for a square pulse 
with zero risetime. In FigfJ (a) we show the results of a ID Boltzmann simulations for a 
square pulse with / = 10 16 W cm -2 , A = 0.25 /im, and a target with n e /n c = 12.5. The 
current density is j y ~ 10 22 c.g.s. units and extends over a distance comparable to d p ~ 
1.2 x 10~ 2 /im. The maximum magnetic field is consistent with Ampere's law, which gives 
B st ~ &'Kjyd p /c ~ 5 MG. Assuming a density n e ~ n Q = 2.2 x 10 23 cm~ 3 for the electrons 
which are instantaneously ionized, one gets a steady velocity v st — j y /en e ~ 10 8 cm s _1 . 
This value is lower than the ejection velocity for hydrogen vi ~ 2 x 10 8 cm s _1 . This suggests 
that effects such as screening, nonzero ionization time, and velocity statistics act to keep the 
steady current well below the values that one may estimate according to the SMM, eq.([TJ). 

Both laser and target parameters where varied in simulations in order to investigate the 
scaling of the magnetic field with them. As an example, Fig.[5] (b) shows the results of a 
simulation for a target of hydrogenic ions with density and thickness identical to Fig.||| (a), 
but where we assumed a nuclear charge Z = 2 and scaled the atomic parameters accordingly 
to x -> Zx, t -> Z 2 t, oj — > Z 2 uj, E — > Z 3 E. In order to have the ionization threshold to 
be exceeded at the same instant, the laser pulse had the same envelope and frequency but 
the intensity was scaled by Z 6 . With respect to the Z = 1 case, we obtain a steady field 
with lower peak amplitude which assumes both positive and negative values. 

We also performed 2D Boltzmann simulations for a pulse obliquely incident at 15° on 
the target. The preliminary results show that the magnetic field is much lower in this case. 
Therefore it appears that the steady magnetic field is sensitive to the interaction geometry. 
In any case, the oblique incidence results further confirm the conclusion that no magnetic 
field capable to affect the transmission through the target is generated. 
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III. OPTICAL MICROCAVITIES AND IONIZATION FRONTS 

A. The model 

In this section, we study effects related to two beam-interference in one spatial dimension 
and for wavelengths in the infrared and optical range. In our numerical experiment, a one- 
dimensional interference pattern is generated via an appropriate "target manufacturing": 
the idea is to place a reflecting mirror on the rear side of the target, the one opposite 
to the laser. Such a mirror might be easily produced by a metallic coating on a glass or 
plastic target. Taking a laser pulse with peak intensity between It/2 and It, being It 
the "threshold" value for ionization, a plasma is produced in the target bulk around the 
maxima of interference pattern produced by the incident wave and the wave reflected at the 
rear mirror. 

Since in this regime we deal only with moderate laser intensities, we may use a simple 
one-dimensional fluid model based on continuity, current and wave equations for an ionizing 
medium, originally proposed by Brunei (1990), modified by the inclusion of the polarization 
current. More details about the model and its validity can be found in Cornolti et al. (1998) 
and Conejero Jarque et al. (1999). 

B. Generation of layered plasmas 

We first consider a target with thickness L = 2ir\, being A = 0.8 /im, and density 
n Q = 10n c . The laser pulse has a sin 2 -shaped envelope with a duration of 80 fs (30 cycles) 
and a peak intensity / = 1.8 x 10 14 W cm -2 . The target parameters are chosen to simulate 
a thin foil solid slab and it is enough to take the density as low as 10n c since the maximum 
electron density always remains much lower than this value. 

The electron density vs. space and time is shown in Fig.f|. A clear layered density pattern 
with a spatial periodicity close to A/2 is produced along nearly all the slab. The layers of 
overdense plasma are produced near the maxima of the interference pattern. These maxima 
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appear at close times because of the effect of the smooth envelope of the laser pulse. The 
resulting quasi-periodic structure of the refractive index has in principle some similarities 
with the widely studied semiconductor microcavities and photonic band-gap materials (see 
reviews by Burstein and Weisbuch (1993) and by Skolnick et al. (1998)). 

C. Optical microcavities 

Since the density in the plasma layers is overcritical, and the layers are created in a 
time shorter than a laser halfcycle, the portions of the standing wave between adjacent 
intensity maxima may be "trapped" into the cavity formed by the two neighboring layers. 
This trapping effect is best seen in the case of a CO2 pulse impinging over a gas target 
with L = A = 10.6 /im and n Q = 5n c ~ 5 x 10 19 cm~ 3 . For this target, two plasma layers 
are produced around the positions x = 0.25A,x = 0.75A. Fig. [5] shows the map of the 
electric field at early (a) and later (b) times, showing the generation of the constructive 
interference pattern which yields the layered ionization (a), and the subsequent trapping of 
the field which remains in the cavity at times longer than the incident pulse duration (b). 
The non-ionized regions between density layers clearly act as optical microcavities. 

Since the microcavity length is L c < A/2, light must have an upshifted wavevector k' > k 
in order to persist inside the cavity. This implies also upshift of the laser frequency with 
oj'l > ujl as seen in Fig.|(b). The upshift decreases the critical density value for the trapped 
radiation and therefore wavelengths much shorter than A escape from the cavity. Due to 
the small fraction of light that tunnels out of the cavity one observes radiation emission 
from the target for a time much longer than the pulse duration. Both the frequency upshift 
and the pulse lengthening may provide experimental diagnostics for microcavity generation. 
The lifetime of the cavities is ultimately limited by processes such as recombination, which 
however should appear on times much longer than the pulse duration of a few tens of 
femtoseconds that are considered here and are available in the laboratory. 
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D. Ionization fronts 

As already shown, in our model target ionization is produced around the maxima of 
the "standing" wave which is generated due to the reflection at the rear mirror. However, 
since ionization is instantaneous on the laser period timescale, it is produced as soon as 
the wave reflected at the rear mirror travels backwards and builds up the standing wave 
by interference. Therefore, a backward propagating ionization front is generated, as seen in 
FigJJ. The density at the front exceeds the critical density. This feature is not obtained for 
a single pulse impinging on a dense target, since it undergoes immediate self- reflection and 
penetrates only in the "skin" surface layer (Macchi et al. 1999). 

An example of "overdense" ionization front is obtained in the case of a CO2 square 
pulse 15 cycles long impinging over a target with n e = 4n c . The n e (x,t) contour plot is 
shown in Fig|| The ionized layers merge into a more homogeneous distribution and a 
"continuous" ionization front appears. The merging appears because the time- and space- 
modulated refractive index perturbs the reflected wave substantially, leading to broadening 
of interference maxima. The velocity of the front in Fig.|6| is near to, or even exceeds at 
some times that of light. This is clearly not a physical "moving mirror" with a velocity 
greater than c, but a reflective surface which is created apparently with such velocity due 
to a space-time phase effect. 

ACKNOWLEDGMENTS 

We acknowledge the scientific contributions of D. Bauer and L. Plaja as well as their 
suggestions. Discussions with G. La Rocca, R. Colombelli, L. Roso, and V. Malyshev are 
also greatly acknowledged. This work has been supported by the European Commission 
through the TMR networks SILASI, contract No. ERBFMRX-CT96-0043, and GAUSEX, 
contract. No. ERBFMRX-CT96-0080. E.C.J, also acknowledges support from the Junta de 
Castilla y Leon (under grant SA56/99). 



S 



REFERENCES 



BAUER, D. 1997 Phys. Rev. A 55, 2180. 

BAUER, D. & MULSER, P. 1999 Phys. Rev. A 59, 569. 

BRUNEL, F. 1990 J. Opt. Soc. Am. B 7, 521. 

BURSTEIN, E. & WEISBUCH, C, eds. 1993 Confined Electrons and Photons. New Physics 
and Applications (NATO ASI Series B: Physics, vol.340, Plenum Press, New York, 1993). 
CONEJERO JARQUE, E., CORNOLTI, F. & MACCHI, A. 2000 J. Phys. B: At. Mol. and 
Opt. Phys. 33, I. 

CORNOLTI, F., MACCHI, A. & CONEJERO JARQUE, E. 1998 in Superstrong Fields in 
Plasmas, Proceedings of the First International Conference (Varenna, Italy, 1997), edited 
by M. Lontano et ai, AIP Conf. Proc. No. 426 (AIP, New York, 1998), p.55. 
GIULIETTI, D., GIZZI, L.A., GIULIETTI, A., MACCHI, A., TEYCHENNE, D., CHESSA, 
P., ROUSSE, A., CHERIAUX, C, CHAMBARET, J.P & DARPENTIGNY, G. 1997 Phys. 
Rev. Lett. 79, 3194. 

MACCHI, A., CONEJERO JARQUE, E., BAUER, D., CORNOLTI, F. & PLAJA, L. 1999 
Phys. Rev. E 59, R36. 

MULSER, P., CORNOLTI, F. & BAUER, D. 1998 Phys. of Plasmas 5, 4466. 
RAE, S. C. & BURNETT, K. 1992 Phys. Rev. A 46, 1084. 

SKOLNICK, M. S., FISHER, T. A. & WHITTAKER D. M. 1998 Semicond. Sci. Technol. 
13, 645. 

TEYCHENNE, D., GIULIETTI, D., GIULIETTI, A. & GIZZI, L. A. 1998 Phys. Rev. E58, 
R1245. 

WILKS, S. C, DAWSON, J. M. & MORI, W. B. 1988 Phys. Rev. Lett. 61, 337. 



9 



FIGURES 
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FIG. 1. Spatial distribution of magnetic field (left) and electron charge density (right) five 
cycles after the end of the pulse, for "sin 2 " pulses of 0.1 a.u. (dotted line), 1 a.u. (dashed line), 
10 a.u. (solid line) maximum amplitude and a "square" pulse of 1 a.u. amplitude (dashed-dotted 
line). All the pulses are 5 cycles long. The electric field atomic unit is E au = 5.1 x 10 9 V cm -1 
(corresponding to / = 3.5 x 10 16 W cm' 2 ). 
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FIG. 2. Grayscale contours of the magnetic field B z (a) and the current density j y (b) five laser 
cycles after the laser pulse end, for a 2D2V Boltzmann simulation. The dashed line in (a) and (b) 
give B z /B a and j y /j , respectively along x = 2\xm. The parameters B Q = 27.7 MG, j = 2.2 x 10 22 
c.g.s. units. The solid lines give neutral density contours. The dashed-dotted lines mark the critical 
density surface. Simulation parameters I = 10 16 W cm~ 2 , A = 0.25/Lim, n e /n c = 12.5. 
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FIG. 3. Profiles of the steady magnetic field B z (solid) and the current j y (dashed) in ID Boltz- 
mann simulations. The parameters common to (a) and (b) are A = 0.25 /xm and n e /n c = 12.5. 
In the case (a) the atomic parameters are those of an hydrogenlike atom with Z = 1, and 
I = 10 16 W cm" 2 , B = 27.7 MG, j Q = 2.2x 10 22 c.g.s. units. In the case (b) Z = 2 and laser param- 
eters are scaled accordingly to x — > Zx, t — > Z 2 t, oj — > Z~ 2 uj, E — > Z~ 3 E; / = 6.4 x 10 17 W cm -2 , 
5 G = 50.6 MG, j = 7 x 10 22 c.g.s. units. 
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FIG. 4. Grayscale contourplot of free electron density n e (x,t) for a "solid" hydrogen target 
with a reflecting "metal" layer on the rear face (thick solid line). The pulse parameters are 
/ = 1.8 x 10 14 W cm" 2 , A = 0.8^m, At L = 30(2vr/o; L ) ~ 80 fs ("sin 2 " envelope). The target 
parameters are L = 2ir\, n Q = 10n c = 1.1 x 10 22 cm~ 3 . 
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FIG. 5. Evolution of the electric field inside the plasma slab during the interaction with the 
incident pulse (a) and 80 cycles later (b). The pulse parameters are I = 1.8 x 10 14 W cm -2 , 
A = 10.6/xm, AtL = 15(2-7!" /Vl) ~ 530 fs ("sin 2 " envelope). The target parameters are L = A, 
n Q = 5n c = 5 x 10 19 cm' 3 . 
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FIG. 6. Grayscale contourplot of n e (x,t) for a hydrogen "gaseous" target with a reflect- 
ing "metal" layer on the right boundary. The pulse has square envelope and parameters 
I = 1.8 x 10 14 W cm' 2 , A = 10.6 urn, At L = 15(2vr/w L ) ~ 530 fs Target parameters are L = 2vrA, 
n e = 4n c = 4 x 10 19 cm' 3 . 
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